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We report the first measurement of inclusive antiproton production at midrapidity in Au + Au col- 
lisions at -Jsnn = 130 GeV by the STAR experiment at RHIC. The antiproton transverse mass dis- 
tributions in the measured transverse momentum range of 0.25 < p± < 0.95 GeV/c are found to fall 
less steeply for more central collisions. The extrapolated antiproton rapidity density is found to scale 
approximately with the negative hadron multiplicity density. 
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We report the first measurement of inclusive anti- 
proton production at midrapidity in Au + Au col- 
lisions at nucleon-nucleon center-of-mass energy of 
■Jsnn = 130 GeV. The measurement was motivated by 
the following: 

(1) Lattice QCD calculations predict that at sufficiently 
high energy density matter should be in a state of de- 
confined quarks and gluons [1]. Large energy densities 
are expected to give rise to increased production of anti- 
baryons relative to lighter mass particles. For example, a 
higher temperature in a (locally) equilibrated system would 
result in a larger relative abundance of antibaryons over 
pions. At high pion density, multiple-pion fusion into 
baryon-antibaryon pairs may contribute significantly to 
the antibaryon yield [2]. Therefore, a measurement of the 
antiproton yield relative to negatively charged hadrons may 
provide information about the energy density reached in 
heavy ion collisions. On the other hand, it has been sug- 
gested in the context of the Skyrme model that the baryon- 
antibaryon production rate can be far above a chemical 
equilibrium estimate [3]. The inclusive antiproton mea- 
surement reported here constitutes an important step to- 
ward the goal of understanding the physics of baryon 
production in heavy ion collisions. 

(2) The mechanism of baryon transport over large ra- 
pidities has been the focus of theoretical investigations 
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[4,5]. The recent measurement of midrapidity antiproton 
to proton (jJ/p) ratio of 0.6 in central Au + Au collisions 
at RHIC [6] indicates that a finite net-baryon number is 
present at midrapidity. This implies that a finite baryon 
number has been transported over 5 units of rapidity in 
these collisions. Transport of incoming baryon number 
over several units of rapidity likely occurs very early in the 
collision and affects the subsequent evolution of the sys- 
tem [4,5,7]. We extract the net-proton multiplicity density 
at midrapidity from the inclusive antiproton measurement 
and the published ~p/ p ratio [6]. 

The measurement reported here was carried out in the 
summer of 2000 at the Relativistic Heavy Ion Collider 
(RHIC) by the STAR (Solenoidal Tracker at RHIC) ex- 
periment. The STAR detector [8] consists of several de- 
tector subsystems in a large solenoid magnet, including a 
time projection chamber (TPC), a scintillator barrel (CTB), 
and two zero degree calorimeters (ZDC) [9]. The mag- 
net was operated at 0.25 T. The CTB measured the en- 
ergy deposited by midrapidity charged particles, and the 
ZDCs measured beamlike neutrons. The coincidence of 
the ZDCs formed the experimental minimum bias trigger, 
and, with an addition of a high CTB signal, provided the 
central collision trigger. 

In the off-line analysis, the collision centrality was de- 
termined from the measured charged particle multiplicity 
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in the pseudorapidity range I17I < 0.75 in the TPC. The 
multiplicity distribution was subdivided into eight central- 
ity bins [10]. The corresponding fractions of events from 
the measured minimum bias data sample are tabulated in 
Table I (first column). The measured minimum bias data 
sample represented about 80% -90% [10] of the theoreti- 
cal total hadronic cross section of 7.2 b [11]. For the two 
most central bins, the central collision trigger events were 
also included. 

Tracks were reconstructed from three-dimensional hits 
in the TPC. The primary interaction point (primary vertex) 
was reconstructed from the tracks. Events with a primary 
vertex within ±30 cm longitudinally of the TPC center 
were used in this analysis. Tracks were required to point 
within 3 cm of the primary vertex (distance of closest ap- 
proach (DCA) cut) and to have at least 25 (of 45 maxi- 
mum possible) hits. The antiproton rapidity was limited 
to |y| < 0.1. Reconstructed momentum was corrected for 
particle energy loss in the detector. Momentum resolution 
for antiprotons was estimated to be about 2% at a trans- 
verse momentum of 0.5 GeV/c. The effect of smearing 
due to the finite momentum resolution on antiproton spec- 
tra is negligible. 

Particle identification was achieved by the measurement 
of the truncated mean energy loss, {dE/dx), of charged 
particles in the TPC gas. At a momentum of 0.5 GeV/c, 
the width of the {dE/dx) distribution for antiprotons was 
found to be about 11% for this analysis. We constructed 
a variable [12], z = ln[{dE/dx)/(dE/dx)w], where 
(dE/dx)w is a single parameter approximation to the 
expected Bethe-Bloch value for antiprotons. Figure 1 
(left panel) shows z versus transverse momentum p L for 
negatively charged particles at midrapidity, demonstrating 
the particle identification capability. We present antipro- 
ton results for p L up to 0.95 GeV/c. Figure 1 (right 
panel) shows the z distributions for two p j_ bins. The z 
distributions were fitted by the sum of three Gaussians 
(nine free parameters) corresponding to different particle 
species. The antiproton raw yield was extracted from the 
fit results for each centrality and p± bin. 

A correction factor was applied to the raw yield to ac- 
count for losses due to acceptance, tracking inefficiency, 
and antiproton absorption in the detector. The overall 




0.5 1 1.5 
P ± (GeV/c) z 

FIG. 1. Left: The z variable versus p ± for negatively charged 
particles. Right: The z distributions for two p± bins, as indi- 
cated by the dashed lines in the left panel. 



reconstruction efficiency including all these effects was 
obtained from a full Monte Carlo (MC) simulation, embed- 
ding MC tracks into real events on the raw data level and 
comparing the reconstructed tracks with the input through 
hit matching in proximity. The overall reconstruction ef- 
ficiency for the most central collisions is about 80% inde- 
pendent of p^ in the range of 0.4 < p± < 0.95 GeV/c 
and about 70% at p L = 0.25 GeV/c. The efficiency in- 
creases with decreasing event multiplicity, by about 10% 
to the most peripheral collisions. 

The antiproton yields reported here include secondary 
products of weak decays. MC studies show that within the 
DCA cut of 3 cm, the overall reconstruction efficiencies as 
a function of the measured p L are identical for secondary 
and primary antiprotons, and the measured antiproton p L 
spectrum is the sum of the primary antiproton p j_ spectrum 
and (0.99 ± 0.05) times the p j_ spectrum of secondary 
antiprotons from weak decays. Secondary antiprotons 
typically carry most of the parent antihyperon p j_ due to 
the decay kinematics, and thus the inclusive antiproton p± 
distribution is similar to the primary antiproton distribu- 
tion when the primary antiproton and antihyperon spectra 
are the same. 

Figure 2 shows the antiproton invariant yield per 
at midrapidity (|y| < 0.1) as a func- 
Here, m± is the transverse mass, 



event, 



2iTm ± dm ± dy ' 

tion of m x — m 



TABLE I. Antiproton fit parameters and yields. Antiprotons are measured at midrapidity (|y| < 0.1) and within 0.25 < p± < 
0.95 GeV/c. Listed errors are statistical. Systematic errors are 7% on the negative hadron dN h -/dr], 10% on the antiproton fiducial 
dN/dy, 15%— 25% on the antiproton total dN/dy estimated from the /? ± -Gaussian parametrization, and 10% on cr p± , T m± , and T s . 



Centrality 




Fiducial dN/dy 


Total dN/dy 




T 




bin 


dN h -/dr] 


(0.25 < p ± < 0.95 GeV/c) 


(pi -Gaussian) 


(MeV) 


(MeV) 


(MeV) 


58%-85% 


17.9 


0.83 ± 0.02 


1.15 ± 0.03 


510 ± 13 


228 ± 11 


190 ± 8 


45%-58% 


47.3 


1.92 ± 0.04 


2.92 ± 0.09 


579 ± 17 


293 ± 17 


233 ± 11 


34% -45% 


78.9 


3.00 ± 0.05 


5.05 ± 0.18 


642 ± 20 


362 ± 22 


274 ± 13 


26%-34% 


115 


4.04 ± 0.07 


7.05 ± 0.28 


668 ± 22 


390 ± 25 


290 ± 14 


18%-26% 


154 


5.29 ± 0.08 


10.44 ± 0.50 


741 ± 27 


481 ± 34 


338 ± 17 


11%-18% 


196 


6.41 ± 0.11 


13.38 ± 0.79 


777 ± 33 


525 ± 43 


359 ± 21 


6%-ll% 


236 


7.84 ± 0.06 


17.19 ± 0.52 


804 ± 17 


568 ± 24 


378 ± 11 


0%-6% 


290 


9.49 ± 0.06 


20.53 ± 0.50 


799 ± 14 


560 ± 19 


374 ± 9 
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FIG. 2. Transverse mass distributions of inclusive antiproton 
invariant yield at midrapidity (\y\ < 0.1). For clarity only five 
centrality bins are shown. Errors shown are statistical. System- 
atic errors are 8% point to point and 10% in the overall normal- 
ization. The solid lines are /? ± -Gaussian fits to the distributions. 



V2 2 
p± + m , and mo is the antiproton mass. The 

uncorrelated point-to-point systematic errors on the spec- 
tra were estimated to be 8% by varying the track cuts and 
by comparing different analysis techniques. Systematic 
error on the overall normalization was estimated to be 
less than 10% by varying the event selection as well as 
the track cuts. As seen in Fig. 2, the spectra in central 
collisions fall less steeply than in peripheral collisions in 
the measured p ± range of 0.25 < p± < 0.95 GeV/c. 

In order to characterize the shape of the spectra 
quantitatively, we fit the spectra to three functional 
forms: (I) Gaussian function in p± [« exp( — p±/2a 2 ± )], 

(II) exponential function in m± [<* exp(— m±/T m± )~\, and 

(III) Boltzmann function [a m± exp(— m±/T s )]. Both 
II and III have been commonly used to characterize m± 
distributions. The Gaussian function in p± can result 
from the Schwinger tunneling mechanism for particle 
production [13,14], but is less commonly used. The three 
functional forms fit the spectra equally well, with similar 
X 2 per degree of freedom of about 1. The p± -Gaussian 
fits are superimposed in Fig. 2. The fit values for cr p± , 
T m± , and T%, reflecting the local spectra shape, are listed 
in Table I as a function of the midrapidity negative hadron 
pseudorapidity density, dN^-fdr] [15]. Systematic errors 
on the fitted values are estimated to be 10% including the 
effect of the point-to-point systematic errors in the spectra. 

Figure 3(a) shows the fitted values for a PL , T m± , and Tb 
as a function of <i/V h - / drj [15], used as a collision central- 




50 100 150 200 250 300 

FIG. 3. Dependence of antiproton fit parameters and yields 
on dN h -/dr). (a) Fitted values of parameters from three 
functional forms described in text: cr PL (I), T m± (II), and r B 
(III), (b) Fiducial yield of inclusive antiprotons in 0.25 < 
p ± < 0.95 GeV/c (points) and the estimated total yield using 
the three functional forms for extrapolation (solid, dashed, 
and dotted), (c) Fiducial yield of net-protons (points) and 
the estimated total net- proton yield using p± -Gaussian for 
extrapolation. Error bars are statistical errors. Systematic errors 
on the fitted parameters in (a) are 10%, and are indicated in 
"caps" in (b) and (c). 

ity estimate. The three parameters clearly exhibit the same 
monotonic trend with centrality: the more central the col- 
lision, the larger the parameters. For comparison, the pro- 
ton and antiproton parameters, T m± , measured in central 
Pb + Pb collisions at the Super Proton Synchrotron (SPS), 
are about 300 MeV [16-18]. This value is similar to our 
result for the peripheral collisions at dN^-fdr] ~ 50. It 
has been suggested that the parameters T m± (and similarly 
a PL and T B ), contain information about transverse radial 
flow which can be generated by a pressure gradient in the 
collision system [16]. The results indicate qualitatively 
that at RHIC we observe a stronger transverse flow for 
midcentral and central events than at SPS. 

We characterize the inclusive antiproton production rate 
by the fiducial rapidity density, dN/dy, in the measured 
p± range, by summing up the data points in each spectrum. 
The results are listed in Table I and shown in Fig. 3(b) as 
a function of dN^-fdr). In order to estimate the total anti- 
proton rapidity density, we extrapolate our measurement 
to all p±. We show in Fig. 3(b) integrals of the three 
fitted functional forms. Comparisons with the negative 
hadron p± distribution [15] reveal that both the m± ex- 
ponential and the Boltzmann extrapolation of our data for 
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central collisions exceed the negative hadron yield around 
p± = 2 GeV/c. Hence, the exponential in m± and the 
Boltzmann distribution are likely to overestimate the anti- 
proton dN/dy. We quote in Table I the p± -Gaussian 
extrapolation as our best estimate of the antiproton total 
dN/dy. By comparing the three integrals, we estimate 
the systematic errors on the quoted total dN/dy to range 
from 5% for the most peripheral bin to 15% for the most 
central bin, in addition to the 10% systematic error on the 
overall normalization. As seen in Fig. 3(b), the total rapid- 
ity density scales approximately with the negative hadron 
multiplicity. However, as the systematic errors are largely 
correlated, there is an indication of a larger ratio of the an- 
tiproton dN/dy to dNh-/drj in central collisions than in 
peripheral collisions. 

The centrality dependence of relative antiproton pro- 
duction is qualitatively different from what has been ob- 
served at lower energies at the SPS (Jsnn — 17 GeV) [18] 
and the alternating gradient synchrotron (^Jsnn ~ 5 GeV) 
[19], where production of antiprotons relative to pions de- 
creases from midcentral to central collisions. Hadronic 
model studies [20,21] show that this decrease is a result 
of a strong absorption of antiprotons in the collision zone, 
but that the initial production of antiprotons relative to pi- 
ons increases with centrality. We note that absorption of 
antiprotons may also play a role at the RHIC energy. 

A constant p~/ p ratio has been measured at midrapidity 
in the p±_ range of 0.4 < p± < 1 GeV/c [6]. Combining 
our results with the p~/ p ratio, we extract the midrapidity 
net-proton (p — ~p) fiducial density within 0.25 < p ± < 
0.95 GeV/c, as shown in Fig. 3(c). The fiducial net- 
proton density increases approximately linearly with the 
negative hadron multiplicity. Figure 3(c) indicates that in 
the most central collisions there are 6.2 ±1.8 protons per 
unit rapidity in excess of antiprotons in the measured p j_ 
range at midrapidity. The total net-proton density at midra- 
pidity can be calculated from the extrapolated antiproton 
yield and the assumption that p/ p is constant over all p L . 
This is represented as the curve in Fig. 3(c). The system- 
atic uncertainties on the total net-proton density are 35%, 
and are largely correlated among centrality bins. There- 
fore, in central collisions approximately 14 net protons per 
unit of rapidity are found at midrapidity. For comparison, 
the HIJING (heavy ion jet interaction generator) model 
predicts a net-proton density of 6 for central Au + Au 
collisions at ^Jsnn = 200 GeV, while the HIJING model 
with the baryon junction mechanism predicts a net-proton 
density of 16 [5,22]. As also seen in Fig. 3(c), the net- 
proton density exhibits a stronger than linear increase with 
the negative hadron multiplicity. The nonlinear increase 
is borne out by the systematic drop of the p~/ p ratio with 
centrality [6] and the slight increase of the antiproton to- 
tal rapidity density. The results indicate that more incom- 
ing baryons, relative to produced particle multiplicity, are 
shifted from beam rapidity to midrapidity for more central 
collisions. 



The inclusive antiproton yield reported here is the sum 
of the primordial antiproton yield and the weak-decay con- 
tributions: p + 0.64(A + 2° + H + n + ) + 0.522". 
We estimate the total antibaryon rapidity density, under 
the assumption of isospin symmetry with n ~ p~ and 
2 = 2 ~ 2 , to be twice the measured antiproton 
rapidity density. Thus, the total net-baryon density is 
approximately twice the total net-proton density presented 
in Fig. 3(c). 

To conclude, we have measured inclusive antiproton 
production at midrapidity (\y\ < 0.1) in the p± range of 
0.25 < p ± < 0.95 GeV/c from 7^ = 130 GeV Au + 
Au collisions at RHIC with the STAR experiment. In the 
measured p± range, the antiproton transverse mass distri- 
butions are found to fall less steeply in more central colli- 
sions. For the most central collisions, the transverse mass 
distribution is significantly flatter than in Pb + Pb colli- 
sions at the SPS. The antiproton rapidity density at midra- 
pidity is found to scale approximately with the negative 
hadron multiplicity density. For the most central collisions, 
over the range 0.25 < p± < 0.95 GeV/c at midrapidity, 
we measure 9.5 ± 1 .0 antiprotons per unit rapidity, result- 
ing in 6.2 ±1.8 protons per unit rapidity in excess of anti- 
protons taking into account our previous measurement of 
the p~/ p ratio. From extrapolation of the yields to all p L , 
we find approximately 14 net protons per unit rapidity at 
midrapidity. 
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